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• Abstract 

, The combined analysis of polarized DIS and SIDIS data is performed in NLO 

CN I QCD. The new parametrization on polarized PDFs is constructed. The uncertain- 

■ ties on PDFs and their first moments are estimated applying the modified Hessian 

. method. The especial attention is paid to the impact of novel SIDIS data on the 

polarized distributions of light sea and strange quarks. In particular, the important 
question of polarized sea symmetry is studied in comparison with the latest results 
on this subject. 

Oh PACS: 13.65.Ni, 13.60.Hb, 13.88.+e 

1^ ! Since the observation of the famous spin crisis in 1987 [1] one of the most intriguing 
and still unsolved problems of the modern high energy physics is the nucleon spin puzzle. 
^ ' The key component of this problem, which attracted the great both theoretical and ex- 
>- ' perimental efforts during many years is the finding of the polarized parton distributions 
^ ■ functions (PDFs) in nucleon. 

fSj ■ The analysis of data on inclusive polarized DIS enables us to extract such important 
^ I quantities as the singlet AS(x, Q^) and nonsinglet Aq3{x,Q^), Aqs{x,Q^) combinations 
00 ■ of the polarized PDFs, and, thereby, the sums of valence and sea PDFs Aq + Aq = 
^ ■ Aqv + 2Aq. Besides, dealing with DIS data, the gluon helicity distribution AG{x, Q"^) 
O ■ is determined due to the evolution in singlet sector and weak dependence on AG of the 
>■ polarized structure function gi in NLO QCD. However, even for singlet combinations 
Am + Am, Ad + Ad, As + As, considered as well-determined within DIS, we still meet 



X 



5-1 the problem: there are only two equations corresponding to inclusive asymmetries Ai 
■ ■ ' measured on proton and deuteron targets from which we should determine three unknown 
combinations Aq^i^x), Aqs{x) and AS(x) (or, alternatively, Au + Au, Ad + Ad, As + 
As). So, it is unavoidable to involve some additional assumptions performing the fitting 
procedure for the purely inclusive DIS data. Moreover, DIS data can not help us to solve 
the important problem of valence and sea PDFs separation. 

The basicQ process which enables us to solve these problems is the process of semi- 
inclusive DIS (SIDIS). However, until recently the quality of the polarized SIDIS data was 
rather poor, so that its inclusion in the analysis did not helped us [2] to solve the main 
task of SIDIS measurements: to extract the polarized sea and valence PDFs of all active 
flavors. Only in 2004 the first polarized SIDIS data with the identification of produced 
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"^As long as no neutrino factory is built and no hyperdense polarized target is created, we cannot study 
the DIS processes with the neutrino beam, which would enable us to find the polarized valence and sea 
quark distributions separately. 
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hadrons (pions and kaons) were published [3]. These data were included in the global 
QCD analysis in Ref. [1]. Further, COMPASS [5] presented the polarized SIDIS data 
(without identification of hadrons in the final state), in particular, in the low x region 
unaccessible for HERMES. This data were included in the latest parametrization of Ref. 
[B]. Recently, the new data on the SIDIS asymmetries A'^ , were published [7] by the 
COMPASS collaboration. It is of importance that this data cover the most important and 
badly investigated low x region. In this paper we include this data in the new global QCD 
analysis of all existing polarized DIS and SIDIS data. The elaborated parametrization on 
the polarized PDFs in some essential points differs from the parametrization of Ref. [6] 
(see below). Nevertheless, the results on PDFs obtained with both parametrizations are 
compatible within the errors. 

Constructing the new parametrization we tried to be as close as possible to our previous 
NLO QCD analysis of pure inclusive DIS data [8]. Namely, to properly describe the DIS 
data we (just as before [8]) parametrize the singlet AS and two nonsinglet Ags, Agg 
combinations at the initial scale Qq = 1 GeV^ 

^ ^ x"^^3 (1 - a;)/^A.3 (1 + ^^^^a;) 



= ^Ags TTZ^ZTT^ ^n/3a,. n , Zm 7E^ ■ j 



^1 - x)^A-?8 (1 + -f^g^X + 6AqsVx) 

■ /(I x"^'8 (1 - x)'5^'8 (1 + 7^^^x + 6AgsV^)dx 
Then, the quantities Au + Au, Ad + Ad, As = As are determined as: 

Au + Au = -(3Ag3 + Agg + 2AE), Ad + Ad = -(SAgg - Agg + 2 AS), (4) 
6 6 

As = As = ^(AS-Agg). (5) 



The gluon PDF is parametrized as 

^ X"^^(1-X)^^^(1 + 7AGX) 

^^^/o^x"AG (1 - xy^o {l + -fAGx)dx ^ ' 

It is easy to see that in our parametrization the coefficients r] are just the first moments of 
the respective local quantities. In particular, the advantage of Ags and Agg parametriza- 
tion in the form ([2]), ([3]) is that with the such choice it is very convenient to apply and 
control the SUf{2) and SUf{3) sum rules 

7] = Aiga = dxAq^{x) = F + D = 1.269 ± 0.003, (7) 

VAgs = Aiqs = C dxAqsix) = 3F - D = 0.586 ± 0.031. (8) 
Jo 

Besides, the such inheritance with Ref. |8] allow us to clearly see the impact of SIDIS 
data on the results of pure inclusive DIS data analysis. 

Further, to properly describe the SIDIS data we, besides AS, Ags and Agg, parametrize 
the sea PDFs of u and d fiavors: 

^OiAu ^'^I^Au _ x"Ad (1 3;)/^Ad 

^''^'^^''jo'x--H^-xr-^dx' ^ "^^^ x-^^1 - x^A^dx 
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Then, Am and Ad are determined from Eqs. (j4]) and ([9]), while the valence PDFs are 
determined by Awy = Am — Am and Ady = Ad — Ad. Thus, all polarized PDFs are 
completely determined within the parametrization. 

Comparing DIS sector, Eqs. ([I])-®, ^ with the respective parametrizations from 
Ref. [H], one can see some distinctions. These are additional factors 7Ag32;, lAgs^ iii Eqs. 

([3]) which are introduced to provide the better flexibility of the parametrizations on 
the respective quantities, required by the inclusion of SIDIS data. Besides, we introduce 
the additional factors 6/\gg a/x in Eq. ([3]) and 'Jagx in Eq. to provide the possibility of 
sign-changing scenarios for As and AG, respectively (see below). We consider two options 
for handling of SUf{2) and SUf{3) sum rules. In the first case we apply the constraints 
([7j) and dHD putting r/Aga and r^Agg equal to the central values of the respective constants. 
In the second case we allow r^Agg and ?7Agg to vary within the uncertainties for F + D and 
3F — D. It is of importance that, as we will see below, both options produce almost the 
same results. 

We analyze the inclusive Ai and semi-inclusive A^ asymmetries. The inclusive asym- 
metry reads 



gip{x,Q'^) 
F2p{x,Q^)/{2x{l + R{x,Q^))}' 



A. = r. ^.wir;. , ... (10) 



where the polarized structure function g^p in NLO QCD looks as 

2g^p = E |Ag + "^^^{AC^ ® Ag + ^AC^ ® AG]| . (11) 



Throughout this paper we use the MS factorization scheme. The respective coefficient 
functions ACq^g can be found in Ref. |9j. The semi- inclusive asymmetries, besides x and 
Q^, depend also on hadronic variable z. As usual, we apply the semi-inclusive asymmetries 

Ah , q2^ ^ Io.2dzgipix,Q^,z) 

^^^'^^ Jl,dzF,%{x,Q\z)/{2x{l + R{x,Q^))} ^ ' 

integrated over the cut z > 0.2, which corresponds to the current fragmentation region. 
The semi-inclusive structure functions g^ in NLO QCD are given by 

K = Ee?Ag[l-h®gAC,,®]D;+(Ee;Ag)®^AC,,®Dj 

<?!<? <?l9 

+ AG®^AC,,®(Ee;/^;). (13) 

q,q 

The respective Wilson coefficients ACgg^gg^gg can be found in [TDj. Of great importance for 
the SIDIS data analysis is the choice of parametrization on the fragmentation functions. 
We use here the latest NLO parametrization from Ref. [11]. Calculating F2 in Eq. 
([inD and F^ in Eq. ([12]) we use parametrization for R from [12] and the recent NLO 
parametrization on unpolarized PDFs from Ref. [TB]. For the as(Q^) calculation we 
apply the same procedure as in Ref. [13] (i.e. as{Q'^) in MS scheme is calculated just 
as in [13] with as^M^) = 0.1145). The deuteron structure functions gid and gi^ are 
calculated applying g['^ = {g[p^ + gin){^ — |'^d)/2 with ujd = 0.058 (see, for instance [S]). 
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In our analysis the positivity constraint |Ag| < q and \^G\ < G holds with the precision 



Of importance is the correct application of the factor 1+7^ = 1 + AM'^x^ /Q"^ in 
the analysis - see the discussion on this question in Ref. pHj. When one rewrites the 
asymmetry Ai in terms of F2 instead of Fi then this factor precisely cancel out in the 
ratio of gi{l + 7^) and Fi = ^2(1 + 7^)/2a;(l + R), so that we just arrive at right-hand 
side of Eq. ffTOl) . The same cancellation holds in SIDIS case, so that one arrives at Eq. 
f fT2|) . At the same time, sometimes [15] the inclusive data is tabulated in the form gi/Fi 
(not in the form Ai = {1 + '~f'^)gi/Fi). In this case we fit the experimental values Qi/Fi 



One of the main conditions of the successful global QCD analysis is the robust program 
for the DGLAP solution. The respective program should be fast, well tested and should 
provide a good precision of DGLAP solution. The elaborated in Ref. [8] program for the 
polarized DIS data analysis, based on inverse Mellin transformation method, satisfies to all 
these requirements. So, we again use here this program package, properly modifying it in 
accordance with the peculiarities of SIDIS data (calculation of SIDIS structure functions 
and asymmetries in the space of Mellin moments are included). 

All procedures of the global QCD analysis are based on the construction of the effective 
function that describes the quality of the fit to data for a given set of varying parameters 
{tti}. In the case of polarized DIS and SIDIS analysis one usually uses the function in 
the form (see, for instance, [6] for detail) 



where Ag^p is the measured value of the asymmetry, 6{Aej.p) is its uncertaintjo, Atheor 
is its theoretical estimation. For the minimization of function we use the MINUIT 
package [16]. 

For our analysis we collected all accessible in literature polarized DIS and SIDIS data. 
We include the inclusive proton data from Refs [TTJ [HI [191 [201 [15] , inclusive deuteron data 
from Refs [2ll [13 [IS [22l [201 [T5j and inclusive neutron data from Refs ^231 [211 [251 [26] . The 
semi-inclusive data (asymmetries A^J^ , A^^^'^ , A^^'^ ) are taken from Refs [3l [271 IS] 
and, besides, we include the latest COMPASS data from Ref [7]. In total we have 232 
points for the inclusive polarized DIS and 202 points for semi-inclusive polarized DIS. If 
we fix ?7Ag3 and rj^g^ by the center values of F + D and 3F — D in sum rules ([7|), ([H|), then 
for 16 fit parameters xllindusive = 188.4 and Xolsemi-indusive = 194.8 for DIS and SIDIS 
data, while Xoltotai = 383.9 for the full set of data (434 points). On the other hand, if t^a^s 
and r/Aqg are varied within the errors on F + D and 3F — D, then for 18 fit parameters the 
resulting values are: xllindusive = 188.2, xl\semi,ndusive = 194.8 with xlltotai = 383.7. 
Thus, one can conclude that the fit quality is quite good: Xo/D.O.F. ~ 0.84. 

The optimal values of our fit parameters are presented in Table [T] for both options, 
with fixed and varied r^Agg and ?7Aqg • As it is seen from Table [H the results are almost the 
same in both cases: the differences in parameters are less than 1%. 

Our calculations show that the fit quality does not decrease if we cancel the extra 
parameters setting /^Agg equal to /^aqs, since their values occur very close to each other 

^We treat {S{Aexp))'^ in Eq. (|14p as the quadratic sum of the statistical and systematical errors (see 
discussion on this question in Ref. |6j). 



0.001. 



by the gi/[{l+r 



2)F2/2x(l + i?)]. 
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Table 1: Optimal values of the global fit parameters at the initial scale Ql = 1 GeV^. 



Parameter 


AS 






a 


1.0227 (1.0216) 


-0.6342 (-0.6380) 


-0.7916 (-0.7827) 


/3 


3.3891 (3.3873) 


3.1418 (3.1398) 




7 


0.0 (fixed) 


23.9180 (24.2240) 


36.8400 (37.1990) 


6 


0.0 (fixed) 


0.0 (fixed) 


-13.7480 (-13.8390) 


7] 


0.3846 (0.3850) 


1.2660 (1.2690) 


0.6170 (0.5860) 


Parameter 


AG 


Am 


Ad 


a 


0.9040 (0.9154) 


-0.3506 (-0.3412) 


0.2802 (0.2852) 


P 


= Pau 


15.0 (fixed) 


= Pau 


7 


-5.6703 (-5.6317) 


0.0 (fixed) 


0.0 (fixed) 


6 


0.0000 (fixed) 


0.0 (fixed) 


0.0 (fixed) 


7] 


-0.1828 (-0.1813 ) 


0.0672 (0.0670) 


-0.0792 (-0.0794) 



when we try to find the best fit. Besides, we use the equality Pau = PaS. = P^g (just 
as in Ref. ^) since the polarized data at x > 0.6 0.7 is just absent, while in the 
region 0.4 -=-0.7 the statistical errors are too large to feel the difference in values of these 
parameters. 

Certainly, the construction of the best fit should be accompanied by the reliable 
method of uncertainties estimation. We choose the modified Hessian method [28], [29] 
which well works (as well as the Lagrange multipliers method - see [6] and references 
therein) even in the case of deviation of profile from the quadratic parabola, and was 
successfully applied in a lot of physical tasks. 

Let us recall that the standard Hessian method is based on the assumption that the 
global is quadratic near the minimum x^: 

M 

Here are the fit parameters values at the global minimum and Hij are the elements of 
the Hessian matrix - the matrix of second derivatives of x^ in the minimum. Then, the 
uncertainty on any physical quantity F which depends on M fit parameters {oj}, can be 
estimated applying 

= Ax^ Y. I^H-^. (16) 

The Hessian method in the simple form f|T6|) is implemented in the MINUIT program 
[T6] ("Hesse procedure") and it perfectly works if x^ has the parabolic profile near the 
minimum. However, in practice we often meet the problems which decrease the reliability 
of this simple version of Hessian method, and main of them is the deviation of x^ from 
parabolic form - see [2H] for detail. To deal with asymmetric x^ profiles we apply the 
modification [28] , [29] of the standard Hessian method, which is proved to be well working 
even in the strongly asymmetric cases. Shortly, the essence of procedure is following. First, 
one starts with still symmetric errors ±5F on the physical quantity F, which one finds via 
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Figure 1: The scheme illustrating the search of asymmetric errors within the modified 
Hessian method. 



Eq. ( |T6| ) rewritten [28] in terms of eigenvectors and eigenvalues of the diagonalized Hessian 
matrixlj. Simultaneously, one calculates the values of fit parameters {oj} corresponding to 
F+S{F) and F — S{F) and, thereby, the respective values. To obtain asymmetric errors 
[29] corresponding to the real profile one varies Ax^ in Eq. (fT6|) (rewritten in terms 
of eigenvectors and eigenvalues) and calculates the respective values of the parameters, 
finding the intersections of x^(-^) curve with the straight line = Xo + ^X^? where Ax^ is 
already fixed quantity (Ax^ = 1 or A^^ = 18.065 here - see below). Then the differences 
of F values in these intersections with F = Fo in the global minimum just determine 
F + 5^^\F) and F — 5^~\F), where 5*^^-*(F) are in general asymmetric uncertainties on 
the physical quantity F - see Fig. [TJ 

Besides, very important question arises about choice of A^^ determining the uncer- 
tainty size. The standard choice is A^^ = 1, just as we did before in Ref. [S] (see also 
[B] and references therein). However, the such choice of Ax^ can lead to underestimation 
of uncertainties, as it was argued in Ref. [31]. The alternative choice of Ax^ (see, for 
example, Refs. [31], [32] and references therein) is based on the equation 

= L 2f(M75)y -p(-5)^^- ('^) 

where P=0.68 (la deviation) is the probability to find the values of all M fitting parame- 
ters inside the hypervolume determined by the condition x^ ^ Xo + ^X^- our case (17 
parameters) the Ax^ value calculated from Eq. f|T7|) is equal to 18.065. We calculate the 
uncertainties for both Ax^ = 1 and Ax^ = 18.065 options. 

The local distributions together with their uncertainties are presented in Fig. [2l The 
first moments of PDFs together with their uncertainties are presented in Table [2J As it 
was mentioned above, the results for both scenarios with fixed and varied rj^g^ and r^Ags 
are almost the same (the differences in fitting parameters are less than 1%). That is why 
we present all figures and Table [2] only for the option with the varied rj/^g^ and r^Agg- 

Let us now discuss the obtained parametrization. First, one can see that the results 
on the first moments AiS = t^ae and AiG = //ag are very close to the respective results 

^For the respective calculations we apply the program ITERATE by J. Pumplin [3D]. 
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x[ A u+Au] X[A d+Ad] 




10-= 10-' 10"' 10"' 10-' 10- 



Figure 2: Best fit PDFs (solid lines) together with their uncertainties for Ax^ = 1 (inner 
bands) and Ax^ = 18.065 (outer bands). 



(scenario with AG < 0) obtained in Ref. [8] in the case of pure inclusive DIS. Indeed, 
even with the more smaller errors, corresponding to the option Ax^ = 1, at = 3 GeV^ 
we have AiS = 0.372tlfil and AiG = -O.lQlt^ofi in this paper and AiS = 0.329lg:[j?^ 
and AiG = —O.ISIIIq^II in Ref. [^, respectively. Notice that the obtained value AiE = 
0.372l|j:E]?^ is even more close to the respective NLO value ao{Q^ = 3GeV^) = 0.35 ±0.03 
obtained by COMPASS [8] directly (without fitting procedure), from the first moment 
of the measured structure function gid- From Table [2] one can see also that in the more 
pessimistic case Ax^ = 18.065, AiG value is just zero within the errors - we will return to 
this question below. The impact of SIDIS data on the AG{x) shape will be also discussed 
some later. What concerns nonsinglet combinations Aq^lx) and Aqs{x), we see (compare 
Table [1] with Table 3 from Ref. [8|) that, even despite the small difference in initial 
scales {Ql = IGeV'^ here and Ql = 3GeV'^ in Ref. [5]), the values of Aq-^^x) are very 
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Table 2: Estimations of the uncertainties on the first moments of polarized PDFs for two 
options of choice. 





Ax' = l 


Ax^ = 18.065 


AS 

Am + Am 
Ad + Ad 
As = As 

AG 

Am 

Ad 


0.3846llj:[i^^^ 

-0.40201^^1 
-0.03871^^^1 

n T OOO+0.0720 

0.06721°:°!?^ 
-0.0792l°:°^i 


0.3846±Hi^ 

o.864o;°:°s^l 

-0.40201^1^^ 
-0.0387381°:°°^^ 

n T OOO+0.1693 

0.0672l°:°?^f 
-0.0792l°:°ig 



similar for both parametrizations. At the same time, the values of Aqs{x) differ much 
more essentially, which occurs due to the impact (first of all on As - see below) of the 
additional SIDIS data. 

Let us now compare our results on AG with the respective results from Refs. [B] 
and [31] . The such comparison seems to be reasonable because among the number of all 
parametrizations applying the DIS data, the SIDIS data is included only in the sequel 
of papers Ref. [2], |6]. At the same time, though SIDIS data is not included in the 
analysis of Ref. [31], the RHIC tt" production data is addecill there (just as in Ref. [B]), 
which should provide significant impact on AG in the RHIC x region [0.05,0.2]. The 
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Figure 3: Obtained parametrization on AG (solid line) in comparison with the respective 
parametrizations from Ref. ^ (dashed line) and Ref. [31j (dot-dashed line) at Qq = 
IGeV^. 

results of comparison are presented in Fig. [31 From this figure one can see that the best 
fit results on AG presented in this paper and in Ref. [6] are very close to each other while 
they are significantly differ (in both AG < and AG > regions) from the respective 
result of Ref. [31], even for the same sign-changing scenario for AG. On the other hand, 
comparing the results on the first moment AiG = 77^, one can see that the central values 
of this quantity are almost the same for all three parametrizations: these are —0.183, 
—0.118 and —0.120 in this paper, Refs. [6] and [31], respectively. Moreover, including 
the uncertainties in comparison (see Table [2] here. Table III in Ref. [B] and Table IV in 

^In this paper we are mainly interested in the study of the SIDIS data impact on the polarized PDFs 
(especially on the still poorly known sea quark PDFs). The influence of RHIC tt" production data on 
our parametrization will be considered some later. 
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XA S=XAS 



0.06^ 
0.04 7 
0.02 7 




-0.02 7 - 
-0.04 7 

-0.06E ^ ^ 

10"^ 10"^ 10"^ 1 

X 

Figure 4: Obtained parametrization on As (solid line) in comparison with the respective 
parametrizations from Ref. [6] (dashed line) and [8] (dot-dashed line), = IGeV"^. 

Ref. [31]), we see that the first moment AiG is just zero within the errors for all three 
parametrizations. Thus, we can conclude that with the present quality of the data it is 
hardly possible to realize is the first moment AiG zero or not. To answer this question 
and also to distinguish between the different shapes of AG{x), still mor^ precise data is 
necessary. 

Let us now discuss the impact of SIDIS data on the polarized strangeness in nucleon. 
It is of importance because, after the appearance of the first results on As extraction 
from SIDIS data performed by HERMES [3], we met the puzzle with the positive As in 
the middle x HERMES region 0.023 < a; < 0.6, while the total moment Ais definitely 
should be negative [33] in accordance with the sum rule ([8]). Thus, to satisfy this require- 
ment, As(x) should possess the compensating negative behavior in the unaccessible for 
HERMES low X region < x < 0.023, i.e. the sign-changing scenario for As should be 
realized. Looking at Fig. [2] we see that this is indeed the case and we produce the best 
fit namely for the sign-changing As scenario, as well as in Ref. [B]. However, we see also 
the distinction in the As shape, in comparison with Ref. [6] - see Fig. |H Namely, while 
within parametrization [6] As changes the sign one time, within our parametrization As 
changes the sign twice. It seems that this distinction occurs due to the inclusion of the 
latest COMPASS semi-inclusive data [7], which allow to better fixl§ As shape. This is 
illustrative to compare the NLO results on As obtained here with the respective results 
of direct As extraction in LO by COMPASS - see Fig. O We see very similar As be- 
havior in both cases. Certainly, one should be careful comparing LO and NLO results. 
However, the LO and NLO results do not differ too drastically, so that the such compar- 
ison is very useful and allow us to make at least qualitative conclusion about the PDFs 
behavior (shape of distributions). The new COMPASS data on the kaon asymmetries in 
the wide Bjorken x region [0.003, 0.7] will be available in the nearest future (the paper in 

^Notice that within this paper we apply the sign-changing scenario for AG{x) and do not consider 
scenario AG{x) > which was considered in Refs. [H] and [3T]. For this scenario the first moment of AG 
is of opposite sign and larger in absolute value, about +0.3-^ +0.4 instead of —0.1 — 0.2. Nevertheless, 
the such scenario also produce quite acceptable Xo value. The such arbitrariness even in the sign of 
AiG once again tell us that we still need more data to properly fix AG and that even inclusion of tt*^ 
production data (Refs. [6], [31]) still poorly enables us to solve the problem. 

^ Notice that after inclusion of the latest COMPASS data to our fit / D.O.F. value becomes small 
only if we allow As to change the sign twice (due to the additional parameter SAqg \/x ) . 
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Figure 5: Obtained NLO parametrization on As (solid line) in comparison with the 
COMPASS results [7j obtained in LO QCD (points with error bars), Q"^ = 3GeV'^. 



preparation). Thus, one can hope to eventually fix the strangeness in nucleon. 

As it was discussed above, the exclusive task for SIDIS is the finding of sea Aq (and, 
consequently, valence Aqv = Aq — Aq) PDFs in nucleon. We again compare our results 
on Au and Ad with the respective results from Ref. [6]. Comparing the central values 
(best fit values) of Au and Ad one can see that they are quite similar - see Fig. [61 At 
the same time, there are also some distinctions, and main of them are connected with the 
sum Au{x) + Ad{x) and its first moment Aiu + Aid. The point is that recently, analyzing 



xAu 



10-' 



10-^ 



10-^ 



xAd 

0.04F^ 
0.03 r 

0.02 r 
0.01 r 


-0.01 r 
-0.02 r 
-0.03 r 
-0.04I 

10 




10 = 



10-^ 



Figure 6: Obtained parametrization on Am and Ad (solid lines) in comparison with the 
respective parametrization of Ref. [6] (dashed lines), = IGeV^. 

the SIDIS data on production COMPASS obtained rather surprising result [5] that 
the sum [Aiu + Aid]{Q'^ = lOGeV^) is just zero within the errors (see Table 2 in Ref. 
0) 

[Ain + Aid]coMPAss = 0.0 ± 0.04 ± 0.03. (18) 

This result was confirmed in the subsequent COMPASS paper [7], where sum Au{x, = 
3 GeV^) + Ad{x, = 3 GeV'^) of the local PDFs was extracted from the measured 
asymmetries Aid,Aj^ ,^id the region 0.004 < x < 0.3 (see Fig. 4 in Ref. [7]) and 
occurs to be about zero in the whole this region (the central values occur in both positive 
and negative vicinities of zero). Thus, at least in the leading order (COMPASS analysis) 
the sum Au{x) + Ad{x) is about zero in the region 3GeV^ < < lOGeV^, which sheds 
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Figure 7: Obtained NLO parametrization on Am + Ad (solid line) in comparison with the 
COMPASS results [7] obtained in LO QCD (points with error bars), = 3 GeV^ 



new light on our understanding of polarized light quark sea. Namely, the sea is extremely 
asymmetric (Am ~ — Aj), on the contrary to the assumption of symmetric sea scenario 
A'u(a:, Qo) = ^dix^Ql)^ applied in the practically alj^ existing parametrizations based 
on the pure inclusive DIS data analysis. Our analysis shows that the sum A-u + Ac? is 
very small quantity in NLO QCD too. It is very illustrative for qualitative comparison to 
put our best NLO QCD fit for Au{x) + Arf(x) (evolved to the COMPASS = ^GeV^) 
on the figure (Fig. 4 in Ref. [?]) showing the COMPASS LO results on this quantity - 
see Fig. [71 From Fig. [7] it is clearly seen that Au{x) + Ad{x) is very small quantity in 
both LO and NLO QCD orders. In turn, the first moment Aim + Aid for the proposed 
parametrization is also just zero within the errors 

[AiM + Ai(J](g2 = iGeV^) = -0.011°:°^, (19) 

even in the case Ax^ = 1. Notice that QCD evolution weakly influences this result. 
Even at extremely large value 100 GeV^ (for instance, for COMPASS upper bound on 
accessible is about 60 GeV^) the quantity Aiu + Aid is still very close to zero within 
the errors, even in the case Ax^ = 1: 

[Aim + Aid] {Q^ = 100 GeV^) = -0.02l°:[!^. (20) 



Notice that looking on parametrization of Ref. [B] (see Fig. 3 there) we see that the 
distributions Au{x) and Ad{x) are also of opposite sign. However, they are significantly 
differ in their absolute values (just as in the DIS parametrization GRSV2000 ^34j - see 
footnote [To]). As a result, the central value of [Aim + Aid]{Ql = IGeV^) in Ref. [6] (see 
Table IV there) is —0.08 (i.e. almost the same as Aid central value —0.11) instead of 
-0.01 in Eq. ([1^]). 

The important remark should be made here. At present, the SIDIS data is of such 
quality that all above conclusions about sea PDFs should be considered as preliminary. 
Indeed, looking at Fig. [2] here and Fig. 3 in Ref. we see that the uncertainties on Am 
and Ad are rather large even for the option Ax^ = 1 , while in the case of Ax^ determined 
by Eq. ([T^ {Ax"^ = 18.065 here) Am and Ad are still the zeros within the errors, and 

The only exception is asymmetric GRSV2000 parametrization [34]. However, though in this 
parametrization Au and Ad are also of opposite sign, they are strongly differ in absolute value. 
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even just to see them within this option for A^^ we need more SIDIS data (first of all the 
expected COMPASS data). 

In conclusion, the new combined analysis of polarized DIS and SIDIS data in NLO 
QCD is presentecJ"]. The impact of modern SIDIS data on polarized PDFs is studied, 
which is of especial importance for the light sea quark PDFs and strangeness in nucleon. 
The obtained results are in agreement with the latest direct leading order COMPASS 
analysis of SIDIS asymmetries [7] as well as with the recent global fit analysis in NLO 
QCD of Ref. [6], where the SIDIS data were also applied. Nevertheless, there also some 
distinctions concerning, first of all, the polarized quark sea peculiarities. At the same time, 
the quality of SIDIS data is still not sufficient to make the eventual conclusions about the 
quantities influenced mainly by SIDIS. In this situation of especial importance becomes 
the direct Aq extraction in NLO QCD, where (just as in LO QCD) the central values of 
asymmetries and their uncertainties directly propagate to the extracted Aq values and 
their errors. The such NLO QCD method, free of any fitting procedure with a lot of 
parameters, was elaborated in the sequel of papers |[33] (see review [35] for details). At 
present the respective analysis of the whole existing polarized DIS and SIDIS data is in 
preparation. In any case, the new COMPASS semi-inclusive data should be available in 
the nearest future. In particular, we expect vr^, data on the proton target in the wide 
X region (today the only available such data is the HERMES data in the narrow x region 
and only for tt^ production), which should essentially increase the precision of Am, Ad 
and As extraction. 

The authors are grateful to Y. Bedfer, E. Kabuss, A. Kotzinian, A. Korzenev, F. 
Kunne, A. Nagaytsev, A. Magnon, M. Stratmann, G. Piragino, E. Rondio, I. Savin, 
H. Santos, A. Sidorov and R. Windmolders for fruitful discussions. The work of O. 
Shevchenko and O. Ivanov was supported by the Russian Foundation for Basic Research 
(Project No. 07-02-01046). 
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